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Abstract 

We discuss the detectability of Warm/Hot Intergalactic medium (WHIM) via the 
absorption lines toward bright point sources with a future X-ray satellite mission, 
XEUS. While we consider bright QSOs as specific examples, the methodology can be 
applied to bright gamma-ray burst afterglows. We create mock absorption spectra 
for bright QSOs (more than 20 QSOs over the all sky) using a light-cone output of a 
cosmological hydrodynamic simulation. We assume that WHIM is under collisional 
and photo- ionization equilibrium. If WHIM has a constant metallicity of Z = 0.1Z Q , 
approximately 2 Ovil absorption line systems with > 3a will be detected on average 
along a random line-of-sight toward bright QSOs up to z = 0.3 for a 30ksec exposure. 
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1. Introduction 

More than 50 percent of cosmic baryons is dark (Fukugita, Hogan, & Peebles 1998; 
Fukugita & Peebles 2004). The abundance, physical properties, and spatial distribution of 
such dark baryons constitute one of the most important unsolved problems in observational 
cosmology. Recent numerical simulations suggest that a major fraction of dark baryons takes a 
form of WHIM (Warm/Hot Intergalactic medium) roughly with temperature 10 5 K < T < 10 7 K 
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and over-density 1 < S < 1000 (Cen & Ostriker 1999a; Dave et al. 2001). WHIM traces the 
large-scale filamentary structure of mass (dark matter) distribution more faithfully than hot 
intracluster gas (T > 10 7 K typically) and galaxies both of which preferentially reside in clusters 
that form around the knot-like intersections of the filamentary regions (e.g., Suto et al. 2004a, 
2004b; Kang et al. 2005). This implies that WHIM carries important cosmological information 
in a complementary fashion to distributions of galaxies in optical and of clusters in X-ray. 

DIOS (Diffuse Intergalactic Oxygen Surveyor) is a dedicated soft X-ray mission which 
aims at unveiling dark baryons in the universe (Ohashi et al. 2004) via oxygen emission from 
WHIM. Detectability of such emission features was examined in detail both for mock blank 
survey (Yoshikawa et al. 2003; hereafter Paper I) and for mock targeted survey of outskirts 
of clusters in the local universe (Yoshikawa et al. 2004; hereafter Paper II). These authors 
find that approximately 20 percent of cosmic baryons leaves emission signatures detectable by 
DIOS. This fraction corresponds to a relatively higher temperature tail (10 6 K < T < 10 7 K) 
of the entire WHIM distribution, and is sensitive to the assumed metallicity of WHIM. They 
adopt a metallicity model proposed by Cen & Ostriker (1999b) and Aguirre et al. (2001), but 
the metallicity of WHIM is admittedly very uncertain at this moment. 

Another complementary strategy of detecting elusive WHIM is to search for its absorp- 
tion signatures. In fact, several independent groups reported successful detections of absorp- 
tion lines of highly ionized metals in spectra of background QSOs and clusters; Ovil and Ovm 
(Fang, Bryan, Canizares 2002), Ovi at z ~ 0.2 in the FUSE spectrum of QSO H1821+643 
(Tripp, Savage Sz Jenkins 2000), and more recently O VII, N VII and Nvi along the line of sight 
toward the Blazar Mkn 421 (Nicastro et al. 2005a; Nicastro et al. 2005b). See also Nicastro 
et al. (2002); Fang et al. (2002); Mathur et al. (2003); Fujimoto et al. (2004). In particular 
Nicastro et al. (2005a) estimate fi^HiM = (2.7+^%) x io-[°/ H l-i adopting a temperature of 
10 61 K for the WHIM, where lO' / 1 ^ 1 is the oxygen-to-hydrogen number density ratio in units 
of 0.1 times the solar value. While the above value seems to be in good agreement with that 
expected from the previous simulations (Cen & Ostriker 1999a; Dave et al. 2001), the estimate 
is still sensitive to the assumed temperature and metallicity of the WHIM along the single line 
of sight to Mkn 421 even apart from its limited statistics. 

This paper examines the detectability of WHIM via absorption, instead of emission, 
lines using the same cosmological hydrodynamic simulation data as Paper I. In particular we 
consider QSOs as a bright X-ray lighthouse, and create mock absorption spectra for a future 
X-ray satellite mission, XEUS (X-ray Evolving-Universe Spectroscopy). This methodology can 
be equally applied to gamma-ray burst (GRB) X-ray afterglows as first proposed by Fiore et 
al. (2000). 

The rest of the paper is organized as follows; section 2 describes our methodology and 
model assumptions in computing mock transmission spectra of the QSOs through WHIM. We 
present examples of mock spectra and their basic statistics in §3. We discuss the detectability 
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of WHIM via absorption in §4 and also consider how to extract physical properties of WHIM. 
Finally section 5 is devoted to the conclusions of the paper and to further discussion. 

2. Creating Mock Transmission Spectra of WHIM toward QSOs from 
Cosmological Hydrodynamic Simulations 

We compute mock transmission spectra of WHIM using the same simulation data as 
Paper I. We create a stacked light-cone output for < z < 0.3 from the simulation data. Then 
we consider 6400 different line-of-sight directions, and compute mock spectra for a bright QSO 
taking into account the cosmic X-ray background (CXB) and UV background (UVB) for photo- 
ionization sources. Our methodology is described in the following subsections. 

2.1. Cosmological Hydrodynamic Simulation 

We use the data of Yoshikawa et al. (2001) who performed a simulation with a hybrid 
code of Particle-Particle-Particle-Mesh (PPPM) Poisson solver and smoothed particle hydro- 
dynamics (SPH). Both dark matter and gas employ 128 3 particles within a periodic simulation 
cube of Lbox = 75/i _1 Mpc per side (in comoving). They adopt a standard ACDM model in 
which the density parameter Q m = 0.3, the baryon density parameter Q b = 0.015/i~ 2 , the di- 
mensionless cosmological constant Q\ = 0.7, the rms density fluctuation smoothed over a scale 
of 8/i~ 1 Mpc, cr 8 = 1.0, and the Hubble constant in units of lOOkm/s/Mpc h = 0.7. The sim- 
ulation includes the effect of radiative cooling, but the energy feedback from supernova and 
the UV background heating are neglected. The simulation that we used here was performed 
in 2000. This is why the adopted value of Qf, = 0.03 is ~ 30% smaller than more recent esti- 
mate, for instance, fib = (0.024 ± 0.001)/i~ 2 (Spergel et al. 2003). In this sense our results may 
underestimate the real absorption signatures by the corresponding amount. 

We create a light-cone output for < z < 0.3 following the procedure described in Paper 
I. The maximum value of z = 0.3 is adopted so as to ensure that the Ovill line is located 
in a spectrum region relatively free from the Galactic confusion (Paper I), while this may be 
very conservative. The output extends a 5° x 5° region on the sky, which corresponds to the 
simulation boxsize at z = 0.3; L box /(l + 0.3) /d\(z = 0.3) ~ 5° with d\(z) being the angular 
diameter distance at redshift z. 

2.2. Calculation of Optical Depths 

Consider one particular line of sight (LOS) with its unit direction vector being n. We 
divide the LOS into 5000 bins Xi {i = 1 ~ 5000) so that their bin width is equal in comoving 
scale, i.e., x i+ i — Xi = dc(z = 0.3)/5000, where dc{z) is the comoving distance at redshift z. The 
gas density p(xj), temperature T(xi) and peculiar velocity v r (xi) of the intergalactic medium 
(IGM) at the 2-th bin are estimated by gathering contribution from all the simulation gas 
particles: 
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N 

p( x i) = ^T j rnW{x i ri-r j ,h j ), (1) 

N rp 

T{x i ) = J2™^W{x i n-r j ,h j ), (2) 
N v r ■ 

v r (xi) = ^2m — W(x i h-r j ,h j ). (3) 
j=i Pi 

In the above expressions, j denotes an index for gas particles ( j = 1 ~ N), Tj , v r j, rj, hj, and 
Pj are temperature, peculiar velocity, position vector, smoothing length, and density of the j-th 
gas particle, and m is the mass of gas particles (independent of j). Yoshikawa et al. (2001) 
classify some fraction of gas particles as "galaxies" if they satisfy both the Jeans condition 
and the over-density criterion. Those particles are excluded from the above summation. We 
adopt the smoothing kernel which is identical to that employed in the simulation (Yoshikawa 
et al. 2001): 

l-(3/2)w 2 + (3/4)w 3 if0<w<l 
= A^\ (2-w) 3 /4 ifl<«<2 , (4) 



otherwise 



where u= |x — v\/h. 

We estimate the number density of an ionized metal A n+ (with charge n+) at the z-th 

bin as 

N /Ha\ X 

n An + ( Xi ) = J2 F A n+ (Tj,pj) I — J — mW{x i YV-r j ,h j ), (5) 

j = i \ n H ' j 

where X is the hydrogen mass fraction (we adopt 0.755), m p is the proton mass, and Z A ,j and 
Fjin+(Tj,pj) represent the metallicity of A and the ionization fraction of A n+ of the j-th gas 
particle. In most cases, we assume a constant metallicity for simplicity: 

Z Atj = 0.1Z A>Q , (6) 

or equivalently a constant number ratio: 

(n A /n H )j = 0.1(n A /n H ) Q (7) 

independently of the time and the local density (see §3.1 below for other models). 

Ionization fraction is calculated using a publicly available routine, SPEX ver. 1.10. 
Papers I and II have assumed collisional ionization equilibrium (CIE) for simplicity because 
detectable emission lines preferentially come from WHIM with relatively high density. In 
the absorption lines, however, lower-density WHIM may become important. So in this paper 
we take account of the photo-ionization effect and compute the ionization fraction of metals 
under the collisional and photo-ionization equilibrium (CPIE) 1 . To do so, we modified the 

1 Non-equilibrium ionization effects on WHIMs have recently considered by Cen & Fang (2006) and Yoshikawa, 
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original SPEX source code, and added routines to calculate the photo- ionization effect. For 
that purpose, we have to specify the differential energy distribution of the diffuse extragalactic 
background radiation. We simply adopt the sum of the cosmic X-ray background and the UV 
background (Miyaji et al. 1998; Shull et al. 1999; Chen et al. 2003): 

J{E) = Jcxb(E) + Jvv(E), (lOeV < E < lOkeV), (8) 

where 

E \- a42 



JcME) = J 1 {^) , (9) 
with Ji = 6.626 x 10~ 26 ergcm -2 s -1 sr -1 Hz _1 (= 1.602 x 10~ 8 ergcm~ 2 s _1 sr _1 keV -1 ), and 

J -( E)=j2 (i3iv)"' 8 - (10) 

with J 2 = 2.4 x 10" 23 ergcm" 2 s" 1 sr^ 1 Hz" 1 (= 5.9 x 10~ 6 ergcm~ 2 s~ 1 sr" 1 keV~ 1 ). At E < 0.26 
keV, Juv is the dominant component of the background radiation, while Jqxb dominates at a 
larger energy region. The resulting ionization fractions of Ovil and Ovill for CIE and CPIE 
are plotted in Figure 1. 

Finally we obtain the optical depth of the ionized metal A n+ along a LOS up to z = 0.3 

as 

rz=0.3 (ft 

r A n+(E) = / n A n+(z)a An +(E') — dz, (11) 
Jz=o dz 

where the rest-frame frequency E' is given by 

E' = E{l + z){l + v r /c). (12) 

Among hundreds of resonance absorption lines for 300eV < E < 2keV (Verner et al. 1996), we 
choose 10 lines for 9 ions that are relevant in the WHIM regime (see Table 1). The absorption 
cross section a A n+ {E 1 ) at E' is computed assuming the Doppler broadening: 



a A n+ {E)=2_, Jk,A n + I — 771 7 6X P 



(E'-E k 



bl El 



(13) 



where k denotes the index of different absorption lines of A n+ , f k) A n + and E k are the oscillator 
strength and the energy of the fc-th line, and br = {2k bT '/m A ) l l 2 is the Doppler fe-parameter 
{m A is the mass of the metal ^4). 

Here we consider the thermal broadening only. Since the natural widths of the lines 
of our interest are less than 0.02eV, we can safely neglect them and thus use Gaussian rather 
than the Voigt profile. Strictly speaking, the internal turbulence is expected to be partly 
incorporated in our simulation data, but may not be accurate enough. Because we are mainly 

Sasaki (2006). The latter paper concluded that the detectability of WHIMs is largely unaffected by the 
non-equilibrium effects, although the line ratio of Ovn and Ovm is systematically overestimated (see their 
Figures 14 and 16). 
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interested in unsaturated lines, the additional effect of the turbulence is not supposed to be 
large. Thus we decided to neglect it in the current analysis. 

2.3. Mock Absorption Spectra of WHIM toward a bright QSO 

Now we are in a position to compute mock absorption spectra of WHIM toward a bright 

QSO. 

Figure 2 shows the cumulative number counts of QSOs in X-ray from ROSAT data 
(Brinkmann, Yuan, Siebert 1997; Yuan et al. 1998). We adopt the flux of /(0.1 - 2.4 keV) = 
7 x 10~ 12 erg s _1 cm -2 as our fiducial value. This value is marked in the vertical dotted line. 
Beyond the value, we find more than 20 QSOs with z > 0.3. For definiteness, we adopt the 
power-law spectrum as 

E \- L5 



F QSO (£) = 1.4 x 10" 12 J erg srW^eV- 1 . (14) 

We consider mock observations of a QSO with 30 ksec exposure. Then fluence T is 
computed by taking account of the absorption due to Galaxy and to the intervening WHIM: 



T(E)=T (E)exp 



J2r An+ (E) 



r-l 



[erg/cm 2 /eV], (15) 



(E)N H ] / E QSO (E)alt [erg/cm 2 /eV], (16) 
Jti 

where Fq(E) determines the observable continuum level free from the WHIM absorption. The 
Galactic absorption is computed using the absorption cross section a ga i(E) of Morrison & 
McCammon (1983) who adopt 1 solar metallicity and the neutral hydrogen column density 
Nh = 3 x 10 20 cm~ 2 . The resulting Galactic extinction factor, exp[— cr ga i(E)NH], is ~ 0.80 at 
E = 0.6keV. For simplicity, we neglect the additional intrinsic absorption. 

2.4- Possible application to bright GRB afterglows 

The same methodology can be applied to bright GRB afterglows. In addition to its 
intrinsic brightness, the more important advantage of GRB for the WHIM science lies in its 
transient nature; once absorption lines due to WHIM are detected, one can perform deep 
follow-up observations to look for the emission counterpart (Papers I and II) after the GRB 
is completely faded. Successful detections of both absorption and emission features will offer 
unambiguous evidence of WHIM, although not so easy in reality. 

Piro (2004) pointed out that the average spectrum energy distribution of a typical GRB 
afterglow is approximated as 

( 1 V s f E Y 



F <"^=Mi000Sj llkevJ ' (17) 
5 =1.2 ±0.2, a =1.13 ±0.07. (18) 

For six years between 1997 and 2002, BeppoSax observed 8% of the entire sky (Jager et al 1997 ). 
The brightest GRB afterglow detected has F = 6 x 10~ 12 erg s~ 1 cm _2 keV _1 approximately 
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(Piro 2004; De Pasquale et al. 2005). 

We consider mock observations of a GRB afterglow for ti <t <tf after the GRB event 
(t = 0). If we observe this GRB afterglow with (ti,tf) = (lday, 1 day + 30ksec), the ffuence is 
comparable with one for our fiducial QSO. For definiteness, we adopt the flux as our fiducial 
value: 

F GM = 6 x 10- (jJL-) ( JL) ^ erg S " WkeV"'. (19) 

Assuming a nominal observational efficiency of 50%, GRB afterglows brighter than our fiducial 
value (Eq.[19]) are roughly expected to be ~ 1/0.5/0.08/6 = 4 per year over the entire sky. In 
reality, however, the above estimate is statistically limited so far, and a better estimate will be 
given from the on-going Swift observation. 

3. Results 

3.1. Identification of absorption lines in mock spectra 

Once mock spectra are constructed, our next task is to attempt the mock observation, 
and in particular, to identify absorption lines in an objective manner. For this purpose, we 
follow the procedure described in Fang, Bryan, Canizares (2002). Consider a schematic example 
of normalized absorption line profiles, exp (— r) as in Figure 3. We set a threshold value 
.Fth = 10 -5 , and identify a contiguous region below the threshold 1 — F t ^ as an absorption line 
system. Suppose the line profile which first down-crosses the 1 — F t h at E — Ei and then up- 
crosses at E = Ej. The equivalent width W (in the observer's frame) of this system centered 
at E = Ei is calculated by integrating 1 — exp (— r) over Ei < E < Ef. 



Wi{Ei) = / l-e- T ^> dE. (20) 

JEi 1 J 

If the line system is identified as a particular species with the oscillator strength the equiva- 
lent width in its unsaturated regime is rewritten in terms of the corresponding column density 
Ni as: 

Ar 9xl0 14 Wi . 

Ni = — cm" 2 . (21) 

1 /, O.leV 1 J V ; 

In the present case, the S/N is determined by the energy resolution of the detector (for 

XEUS, the effective area S c s = 60000cm 2 and the energy resolution AE = 2eV; for comparison, 

intrinsic widths of the lines of our interest are typically below 1 eV). Therefore the S/N at the 

energy of the line center, Ei, is estimated as 

„,., (rJ mE,)-M(E l )]AE 

simply from the Poisson statistics of the observed photons, where M and Mo are the photon 
number counts [ph/eV] with absorption and without absorption (continuum level), respectively. 
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We compute the equivalent width and the S/N of absorption line systems for 6400 
random LOSs in the 5° x 5° region assuming the fluence described in §2.3. Figure 4 plots the 
cumulative number distribution of the equivalent width of absorption line systems per LOS 
(z = 0.0 — 0.3) from our mock observation for 10 strong lines (Table 1). On average, O VII (574 
eV) is the most prominent line in terms of the equivalent width, and O VIII (654 eV) and Fe 
xvn (826 eV) are the next. Figure 5 shows the corresponding cumulative number distribution 
of S/N. Note that the value of S/N should be replaced by 

S/N = (S/N) Mncial ^ 6xlQ _ 12 [erg 8 - W2kE v-i] J 

-(VN) / W = 500eV) \ 1/2 m) 
- (S/iV) fiducial y 8800[ph/eV] ) (23) 

if Fq, and more generally the number of photons during the exposure time, A/"o, are different 
from our fiducial values. In our adopted model for the QSO and the XEUS performance, O VII 
574 eV and Ovm 654 eV lines achieve S/N = 3 when their equivalent widths are W = 0.05 eV. 
We take this value as our detection limit. 

We find that the number of absorption systems which have S/N > 3 is 1.58/LOS (Ovil 
574 eV) and 0.37/LOS (O vm 654 eV) as summarized in Table 2. The other eight lines exhibit 
below 0.1 systems/LOS. The number of the case that both Ovil (574 eV) and Ovm (654 eV) 
exhibit S/N > 3 at the same position of WHIM is 0.35/LOS. 

In order to check the dependence on the assumed metallicity model, we also consider two 
other metallicity models. One is a phenomenological model adopted in Paper II (the radiation 
pressure ejection model of Aguirre et al. 2001): 

Z A>j /Z A>Q = min [0.2,0.02(p,/p b ) a3 ] , (24) 

where pb is the mean baryon density at redshift z. The other adopts a constant metallicity 
comparable to the typical value for intra-cluster medium: 

Z Aj = 0.3Z AQ . (25) 

Figure 6 shows the dependence of the cumulative number distribution of oxygen absorption line 
systems on the metallicity models; O VII (574eV) and O vm (654eV) in thick and thin lines, 
respectively. Note that Z A j > 0.1Z A<Q is satisfied only in relatively dense regions Pj/pb > 214 
(eq.[24]). This is why our fiducial model (Z A j = O.IZ40) has more absorption line systems 
than the model of equation (24). 

We also show the results toward our fiducial GRB afterglow in Table 3. These are almost 
equivalent to those toward our fiducial QSO (Table 2). Thus, we concentrate on the results 
toward QSOs in what follows. 
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3.2. Detectability of Oxygen Absorption Line Systems 

Since the most prominent absorption lines detectable by XEUS are Ovil and Ovm, we 
focus on these lines in what follows. Figure 7 shows column density maps of Ovil {Left panel) 
and Ovm {Right panel) for < z < 0.3. As indicated in those panels, we select three different 
LOSs (A, B, and C). The mock spectra for those LOSs are shown in Figure 8; spectrum A 
exhibits many prominent absorption lines, and spectrum B has one strong line system while 
spectrum C has none. 

We find two distinct absorption line systems which are located at z — 0.12 (marked by 
f) and z = 0.29 (marked by i.) along the LOS of A; Ovil {z = 0.29, E = 445eV, W = 0.16eV), 
Ovm {z = 0.29, E = 507eV,W = 0.09eV), Ovn {z = 0.12, E = 514eV, W = O.lOeV), and Ovm 
{z = 0.12, E = 585eV, W = 0.12eV). Figure 9 plots the column density map of Ovil {Left panels) 
and Ovm {Right panels) but for 0.11 < z < 0.14 ( Upper panels) and for 0.26 < z < 0.30 {Lower 
panels) so as to depict the WHIM responsible for the absorption systems in the spectrum A. 

Figure 10 shows the mass-weighted (within each cell) temperature and the over-density 
distributions along the three LOSs. The top two panels for the region A have high over-density 
(10-100) and high temperature (10 6 - 10 7 K) WHIM clumps at z = 0.12 and z = 0.29. They are 
indeed responsible for the prominent lines in the spectrum A of Figure 8. Spectrum B has one 
absorption system at z = 0.22 (Ovil 574 eV W = 0.06eV) as indicated in Figures 9 and 10. In 
reality, however, the identification of the line redshift is very difficult from a single line alone. 
One may assume statistically that the strongest line in a single LOS should be that of Ovil 
574 eV, and then infer the redshift 574/470 — 1 = 0.22. In this particular example, it turns out 
a correct guess, but may not work always. Finally spectrum C has no prominent absorption 
feature. Indeed Figure 10 shows that there is no region with over-density 5 = p/p > 10 along 
the LOS. 

Figure 11 plots the distributions of mass- weighted temperatures and over-densities of 
WHIM averaged over each absorption line system of Ovil (574 eV) and O VIII (654 eV). 
Spatial boundaries of absorption line systems are somewhat ambiguous. We define here the 
cosmic mean density as the threshold as illustrated in Figure 12, and we identify one connected 
region exceeding 5 = 1 as a single absorption system even if there are many density peaks and 
valleys inside. Then we compute its temperature and over-density by mass-weighted averaging 
over the relevant SPH particles within the entire line system. 

The different symbols indicate the range of the equivalent width W of the absorption 
line systems; circles, crosses and triangles correspond to the systems with W > 0.12 eV, 0.07 < 
W < 0.12 eV, and 0.05 <W< 0.07 eV, respectively. Both Ovil and O VIII can probe WHIM 
with T > 10 6 — 10 7 K and over-density of 5 = 10 to 1000. WHIM with T < 10 6 K, however, can 
be traced only by O VII. 

Figure 13 plots the fraction of gas contained in oxygen absorption line systems (O VII 
in Left panels and O vm in Right panels) which have the equivalent width larger than the 
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specified value W. In practice we first select all oxygen absorption line systems along the i-th 
LOS whose equivalent width exceeds W, and sum up their hydrogen gas column densities. We 
denote the sum as iVf}™(> W). Then we take the following ratio averaged over 6400 LOSs: 

6400 

R h{> W) = ^ (26) 

E^(>°) 
i=i 

for O VII and O VIII separately. In reality this is an unobservable quantity, but provides a 
simulated gas mass fraction of oxygen absorption line systems which have a detectable line 
width. More importantly, this is the essential factor in properly estimating the total cosmic 
baryon density from the detected fraction of oxygens (§4.2). 

In the upper panels, three shaded regions correspond to the temperature ranges T > 10 7 
K, 10 6 < T < 10 7 K, and T < 10 6 K from top to bottom, and in the lower panels, three shaded 
regions correspond to the density contrast ranges 5 > 100, 10 < 5 < 100,5 < 10 from top to 
bottom. Comparison of the left and right panels suggests that Ovill traces higher temperature 
and denser regions than O VII. In the case of O VII (574 eV), we can detect only 32 % of 
baryon above the detection limit (W ~ 0.05 eV; §3.1). The region of T > 10 7 K accounts 
for 19 % of detectable baryons, which may be located in virialized regions like intra-cluster 
medium and thus may be already detected through their thermal bremsstrahlung. The regions 
of 10 6 < T < 10 7 K and T < 10 6 K account for 68 % and 12 % of detectable baryons, which can 
be regarded as WHIM. Note that one can detect the absorption lines of O VII (574eV) with 
T < 10 6 K, while it is almost impossible to find their emission counterparts (Paper I). 

4. Physical Properties of the identified WHIM 

4-1. Ovn and Ovill Absorption Lines 

As we have seen, the most prominent absorption lines are Ovn (574 eV) and Ovill (654 
eV). On average, 35 % of random LOSs have an intervening clump of WHIM which exhibits 
both absorption lines with S/N > 3. In this case, the ratio of their equivalent widths should 
provide the temperature of the WHIM in principle. Consider a clump of WHIM which has a 
uniform temperature T and density p, then the ratio of the equivalent widths for unsaturated 
lines that are relevant here is written as 

H / OVII,574eV /oVII,574eV-Pbvil(^) p) 



(27) 

H / OVIII,654eV /oVIII,654eV-fbviIl(r, p) 

where / is the oscillator strength and F is the ionization fraction. The WHIM which exhibits 
Ovill (654 eV) absorption has temperature of 10 6 — 10 7 K (Figure 11). In this temperature 
range, the ionization fraction approaches the value in CIE, and the ratio of the oxygen equivalent 
widths, eq.(27), is expected to be insensitive to the density, p. 
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Figure 14 provides the scatter plot of the ratio of O VII and O VIII equivalent width 
versus the mass-weighted average temperature for systems where which both VII (574 eV) 
and vill (654 eV) exhibit S/N > 3. Three solid lines indicate the theoretical curves in Eq. 
(27) applying different hydrogen number density nn = 10 -5 , 10 -6 , and 10 -7 cm -3 from top to 
bottom. Inhomogeneous density and temperature structures are responsible for the departure 
from the theoretical model which assumes uniform distribution. Figure 15 shows an example of 
sub-structure in a simulated absorption line system. Temperature is the highest near the center 
of the clump, and over-density is also high. Therefore, ionization fractions of Ovil and Ovm 
are rather suppressed. This anti-correlation of oxygen abundances inside the inhomogeneous 
distribution of WHIMs explains why the simulated systems in Figure 14 preferentially lie at 
higher temperature regime than the theoretical prediction. 

4-2. Estimating from WHIM observations 

If absorption lines of WHIM are detected along a particular line of sight, one would 
naturally attempt to estimate the contribution of WHIM to the total cosmic baryon density, 
qWhim Nicastro et a i. (2005a) indeed reported fi^ HIM = (2.7+^%) x 10-1°/^ using two Ovil 
absorption lines. We perform statistical analysis of f2^ HIM estimation using our mock data. 

Let us first define the estimated fraction of WHIM with signal-to-noise ratio exceeding 
a given threshold, S/N, along one LOS as follows: 

WHIM,est/. C//\n - 1 N H%> S/N) 

n - (> S/N) ~ d c ( z = o.3 )Pc V x(i + Zl y ' (28) 

rest I v. 

0/AA N vii,i{> S/N) fn H 



NZ(> S/N) = ' ' — (29) 

where % denotes an index for absorption line systems, X{= 0.755) is the hydrogen mass fraction, 
and d c (z = 0.3) is the comoving distance from z = to z = 0.3. We convert the observed 
equivalent width into A r ovn(> S/N) using Eq. (21) assuming that the line is unsaturated, i.e., 
in its linear regime of the growth curve. At present, no simultaneous detection of Ovil and Ovm 
absorption lines in WHIM is reported. Consider first, however, the case in which both O VII and 
O Vill are detected. Then we estimate temperature T using Eq. (27) and derive the ionization 
fraction Fqyh(T) assuming collisional ionization equilibrium. If we multiply the hydrogen-to- 
oxygen number ratio, we have an estimate of the hydrogen column density for the z-th absorption 
line system, Nfj\(> S/N). In practice, we adopt (no/n H ) = 0.1(n o /n H ) Q = 10~ 4,07 (Anders & 
Grevesse 1989). Finally this is converted to the density parameter of gas, f2^ IIM ' cst (> S/N) 
responsible for absorption exceeding the detection S/N. 

As we have seen in the previous sections, the inhomogeneity of the temperature and 
density within one absorption system would result in a systematic underestimate of Nff\(> 
S/N). In order to check the systematic effect as well as the assumption of CIE in equation 
(29) (while we did adopt CPIE in creating mock spectra from simulation data), we compare 
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n h%> S/N) with N$(> S/N) directly computed from simulation data (§3.2). Figure 16 
clearly shows the degree of the underestimate bias. 

In order to estimate the total baryon density, fi| st , from r2^ IM ' cst (> S/N), we have to 
correct for (i) the fraction of gas below the detection limit, (ii) the fraction of stars in galaxies, 
and (iii) the above systematic bias. Therefore we write the estimate as 

Q WHIM,est (> S/N) /gas + /star jygm (> g/N) v 

b R H (>S/N) / gas \N^(>S/N)r 1 ) 

In practice, we choose S/N = 3 in the following. We compute the fraction of WHIMs which 
exhibit both Ovil and Ovill absorptions (similarly as in Figure 13), and find Rh(> 3cr) = 0.22. 
As for the second factor, we evaluate directly from simulations; / gas = 0.79 and / sta r = 0.21. As 
Figure 16 indicates, the last correction factor significantly varies from line to line, and it is not 
easy to choose a relevant average value. If we adopt the mean value as plotted in triangles, 
(iVf} m (> 3a)/A^|f t (> 3er)) ~ 1/4. Adopting those correction factors, we plot the histogram of 
Qest - n pjg Ure Y! 

Upper panel is the case where both Ovil and Ovm absorptions are detected (1903 LOSs 
out of 6400 LOSs). The mean value of the resulting f2£ st averaged over the entire LOSs (even 
including the non-detected 4097 LOSs) amounts to 0.023 (plotted in the dotted line). This is 
in reasonable agreement with the value adopted in the simulation, 0.03. On the other hand, 
it clearly demonstrates the fact that any estimate based on the positive detection alone may 
significantly overestimate the real value; it is essential to take account of the above correction 
factors properly. 

Lower panel in Figure 17 indicates the 5217 LOSs along which only O VII absorption lines 
are detected. In this case, one has to assume the temperature a priori, and we show results 
for T = 10 5 ' 5 , 10 6 , and 10 6,5 K. The temperature dependence is very strong and even not 
monotonic. This reflects the strong dependence of Ovil fraction on temperature as exhibited 
in Figure 1. The two peaks in the histogram correspond to LOSs which have single and 
double absorption line systems. If we assume T = 10 6 K (the ionization fraction of Ovil peaks 
around this temperature), Ovil absorption lines alone provides a fairly accurate estimate for 
fi b (Nicastro et al. 2005a). 

5. Conclusions and Discussion 

In this paper, we have examined the detectability of WHIM via the absorption lines in 
a bright QSO spectra using the cosmological hydrodynamic simulation. We created the mock 
XEUS absorption spectra toward a bright QSO assuming the collisional and photo-ionization 
equilibrium. We conclude that on average XEUS will detect 1.58 (O VII 574 eV) and 0.37 
(Ovm 654 eV) per a random line-of-sight up to z = 0.3 (S/N > 3). 

If both Ovil and Ovm absorptions are detected (~ 40% chance) for the same WHIM 
clump, one can estimate its temperature from their line ratios and attempt to infer the cosmic 
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baryon density. As we discussed in detail, however, the reliable estimate requires several careful 
correction factors and still the resulting distribution function of the estimate is fairly wide. So 
any attempt on the basis of small statistical samples should be interpreted with caution. 

In passing we have to emphasize once again that our simulation data adopted Qb = 0.03 
rather than the currently more favored value 0.04. Also we restrict our current analysis up to 
z = 0.3 where a spectrum is relatively free from the Galactic confusion, but in principle one 
can explore WHIM systems at higher redshifts (z ~ 1) with QSO(GRB). Therefore while our 
present conclusions concerning the detectability of WHIM in absorption toward a QSO(GRB) 
may seem fairly modest, most likely they are very conservative. Moreover the result is very 
sensitive to the metallicity of WHIM which is rather uncertain and is likely to significantly 
vary from place to place. Given the strong metallicity dependence as shown in Figure 6, what 
is relevant is not the mean metallicity of WHIM but the fraction of WHIM clumps which 
exhibit detectable equivalent width. In this respect, the wide variation of the metallicity, even 
if its overall mean is O.1Z as we most adopted here, is expected to systematically increase the 
detectability. 

Throughout this paper, we adopted AE = 2eV with calorimeter on-board XEUS in mind. 
Current observations of WHIM by several grating detectors are indicative. With XEUS's large 
effective area, WHIM absorption features will be collected with a much larger number of samples 
as a by-product of QSO studies. With careful statical studies as we presented in this paper, 
these significantly contribute to the proper understanding of physics of WHIM. If higher energy 
resolution with gratings, for example, AE = 0.15eV at E — 0.6keV, becomes available in future, 
it is even possible to obtain line profile accurately and estimate the effect of turbulence. Also 
one can study multiphase nature in a single absorption system. Furthermore, since signal-to- 
noise ratio improves as 1/ V AE, one can significantly increase the number of detectable oxygen 
and other ion absorption systems. These significantly contribute to the proper understanding 
of physics of WHIM (Nicastro 2005). 

In order to be more quantitative, we need much more advanced simulation datasets which 
include more realistic physical processes such as radiative cooling, photo-ionization, star forma- 
tion, supernova feedback, and metallicity evolution. This is exactly what we are performing, 
and we hope to revisit the above issues in due course. 

Numerical computations presented in this paper were carried out at ADAC (the 
Astronomical Data Analysis Center) of the National Astronomical Observatory, Japan (project 
ID: mky05a). We thank the referee Fabrizio Nicastro for very constructive and useful com- 
ments. We believe that his numerous suggestions significantly improve the present paper. We 
are also grateful to Kei Hiraki and Mari Inaba for their generous allocation of computer re- 
sources at the University of Tokyo supported by the Special Coordination Fund for Promoting 
Science and Technology, Ministry of Education, Culture, Sport, Science and Technology. This 
work was supported in part by Grants-in-Aid for Scientific Research from the Japan Society 
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for Promotion of Science (Nos.14102004, 14204017, 15340088, 15740157, and 16340053). 
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Table 1. Ten prominent resonant absorption lines in 0.3 - 2.0 keV 



species 


energy 


oscillator strength 


Cvi 


368 eV 


0.416 


Nvi 


431 eV 


0.675 


Nvn 


500 eV 


0.416 


Ovn 


574 eV 


0.696 


Ovn 


666 eV 


0.146 


Ovin 


654 eV 


0.416 


Neix 


922 eV 


0.724 


Nex 


1022 eV 


0.416 


Mgxi 


1352 eV 


0.742 


Fexvn 


826 eV 


2.96 



Table 2. Expected number of oxygen absorption line systems with S/N > 3 and S/N > 2 per LOS up to z = 0.3 toward a 
QSO discussed in §2.3. 









expected number/LOS 




metallicity model 




Ovn (574 


eV) Ovin (654 eV) Ovn and Ovm 


Z = 0.1Z Q 


S/N>3 


1.58 


0.37 


0.35 




S/N>2 


3.07 


0.82 


0.79 


Z = min[O.2,O.O2(p/p)°- 3 ]Z 


S/N>3 


0.50 


0.15 


0.14 




S/N>2 


0.97 


0.31 


0.28 


Z = 0.3Z© 


S/N>3 


6.71 


2.19 


2.10 




S/N>2 


9.91 


3.40 


3.22 
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Table 3. Expected number of oxygen absorption line systems with S/N > 3 and S/N > 2 per LOS up to z = 0.3 toward a 
GRB afterglow discussed in §2.4. 









expected number/LOS 




mctallicity model 




Ovu (574 


eV) Ovin (654 eV) Ovu and Ovin 


Z = 0.1Z G 


5/iV>3 


1.71 


0.43 


0.41 




S/N>2 


3.30 


0.94 


0.90 


Z = min[O.2,O.O2(p/p)°- 3 ]Z 


S/N>3 


0.55 


0.17 


0.15 




S/N>2 


1.05 


0.35 


0.32 


Z = O.3Z 


S/N>3 


7.08 


2.39 


2.28 




S/N>2 


10.3 


3.66 


3.45 
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T[K] T[K] 

Fig. 1. Ionization fractions of Ovu and Ovin as a function of temperature. Solid 
lines represent results for CIE (Collisional Ionization Equilibrium), while other dif- 
ferent lines indicate those for CPIE (Collisional and Photo-Ionization Equilibrium) 
adopting the hydrogen number density nn = 10~ 4 , 10~ 5 , 10~ 6 , and 1CP 7 cm~ 3 . 
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FLUX (0.1-2.4 keV) [erg s ' cm z ] 

Fig. 2. Cumulative number counts of QSOs in X-ray from ROSAT data as a function of their 0.1 - 2. 
keV flux. Solid and dashed lines indicate QSOs with z > 0.1 and z > 0.3, respectively. 
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0.01 0.02 0.05 0.1 



W [eV] 

Fig. 4. Cumulative distribution of expected number of absorption line systems for < z < 0.3 as a function 
of their equivalent width; Upper panel: Z = O.IZq, Lower panel: Z/Zq = min[0.2, 0.02(p//?b) ' 3 ]- Solid 
vertical lines indicate our fiducial detection limit for the equivalent width (W = 0.05eV). 



21 



o 




S/N 

Fig. 5. Cumulative distribution of expected number of absorption line systems for < z < 0.3 as a function 
of their detection signal-to-noise ratios. The plots are converted from those of Fig. 4 assuming the XEUS 
observation toward a QSO discussed in §2.3. Solid vertical lines indicate our fiducial detection limit 
(S/N = 3). 
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Fig. 6. Cumulative distribution of Ovu and Ovni absorption lines for different metallicity models as a 
function of W (Upper panel) and S/N (Lower panel). Dotted vertical lines indicate our fiducial detection 
limits (W = 0.05eV and S/N = 3). 
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Fig. 7. Column density maps of Ovn (Left) and Ovni (Right) for 0.0 < z < 0.3. Three red squares labeled 
A, B, and C indicate the regions corresponding to Figures 8 to 10. 



24 




300 400 500 600 700 800 

E [eV] 

Fig. 8. Mock transmission spectra of a QSO. Top panel plots fmences through three different LOSs (A, 
B, and C indicated in Figure 7). For an illustrative purpose, curves B and C are artificially multiplied 
by a factor of 10 -0 ' 2 and 10 -0 ' 4 . Lower panels show their residual flucnce normalized by the continuum 
level. The metal line systems with labels f and J in the spectrum A correspond to the WHIM clumps at 
z = 0.12 ( Upper panels in Fig. 9) and z = 0.29 (Lower panels in Fig. 9), respectively. 
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X [deg] X [deg] 

Fig. 9. Column density maps of Ovn (Left) and Ovin (Right) for for 0.11 < z < 0.14 (Upper) and for 
0.26<z<0.30 (Lower). 
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Fig. 10. Rcdshift distribution of temperature and over-density along the three LOSs A, B, and C (from 
top to bottom). 
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Fig. 11. Scatter plots of over-density and temperature of absorption line systems; Ovn 574 eV (Upper 
panel) and O vili 654 eV (Lower panel). Circles, crosses, and triangles correspond to the systems with 
equivalent width of W > 0.12 eV, 0.07 < W < 0.12 eV, and 0.05 < W < 0.07 eV, respectively. 
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Fig. 12. Schematic identification method of clumps of WHIM. 
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Fig. 13. Cumulative ratio of hydrogen column densities (eq.(26)) as a function of the equivalent width 
of the detected absorption line systems; Ovil 574 cV (Left) and O vin 654 eV (Right). The three shaded 
regions indicate the different temperature (Upper) and density (Lower) ranges; T > 10 7 K, 10 6 < T < 10' 
K, and T < 10 5 K ( Upper), and S > 100, 10 < 5 < 100, and 5 < 10 (Lower) from top to bottom. 
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" OVIK " OVI 1 1 

Fig. 14. Ratio of the Ovn (574 cV) and O vin (654 cV) equivalent widths plotted against the 
mass-weighted average temperature. The solid curves show the theoretical curves in the uniform case 
under CPIE assuming the hydrogen number density nn = 10~ 5 ,10 -6 , and 10~ 7 cm~ 3 from top to bottom. 
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Fig. 15. An example of sub-structure in a WHIM which forms a single absorption line system according 
to our identification procedure. From top to bottom plotted are temperature, over-density, number density 
of vn, and number density of VIII. 
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Fig. 16. Comparison between A r |f t (> 3cr) and N^ m (> 3a). Each dot represents a single WHIM clump 
which exhibits both Ovn and Ovni with S/N > 3 (1903 clumps in total). The filled circles and triangles 
indicate the mean and median values of N^ m (> 3a) for the fixed bin with respect to A r ]?f t (> 3a). The 
solid, dotted, and dashed lines show iVjf(> 3a)/N s j n (> 3a) = 1, 1/2, and 1/4, respectively. 
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Fig. 17. Histogram of £!£ st estimated according to eq.(30). The upper panel shows the case that both 
Ovn and Ovni are detected and temperature is estimated from their line ratio. The vertical dotted line 
is the mean value 0.023 computed from the entire LOSs including non-detection. The lower panel shows 
an estimation using only O vn absorption lines assuming temperature T = 10 5 5 , 10 6 , and 10 6 5 K. The 
arrow indicates the value of r2j,(— 0.03) that is adopted in our numerical simulation. 
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